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Abstract 



Using an a decay level scheme, the fine structure in odd nuclei is ex- 
plained by taking into account the radial and rotational couplings between 
the unpaired valence nucleon and the core of the decaying system. It is shown 
that the experimental behavior of the a decay fine structure phenomenon is 
governed by the dynamical characteristics of the system. 
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The a decay fine structure was discovered by Rosenblum j|] since 1929 by measuring the 
ranges of the emitted particle in the air. Usually, theoretical attempts to investigate this 
phenomenon are based on the calculation of the overlaps between, on one hand, the ground- 
state wave function of the parent and, on the other hand, the antisymmetric product between 
the wave functions of the nascent fragments in different configurations after the scission [^,|3| . 
However, quantitatively this phenomenon was not explained rigorously 

It was evidenced ||, at least formally, that the significance of the emitted particle pre- 
formation and that of its penetrability through the barrier calculated from the ground-state 
configuration of the parent up to the scission point are equivalent. This equivalence gives a 
support in the attempt to investigate the a decay process using fission theories. 

Recently, a theory based on the Landau-Zener effect was developed |5|-[7| intending to 
describe quantitatively the cluster decay fine structure phenomenon. This new formalism 
succeeds to reproduce the values of the fine structure hindrance factors for the 14 C emis- 
sion from 223 Ra being in a better agreement than other microscopic theories || with the 
experiment. In this alternative description, the cluster decay fine structure is caused by 
the promotion of the valence unpaired nucleon on some excited daughter levels during the 
disintegration process. It was claimed that the same promotion effect can also govern the 
fine structure in the case of a decay. The first step in such a treatment is the elaboration 
of a two-center realistic level diagram during the whole disintegration process, that means, 
starting from the parent single-particle energy distribution and following the variations of 
the level energies up to their asymptotic configuration attributed to the separated daughter 
nucleus and the alpha particle. In this picture, it is only intended to treat the a cluster 
with a smooth potential in order to estimate its influence upon the levels of the daughter 
during the disintegration process, like a polarization effect. The particle-core couplings are 
produced merely between levels belonging to the nascent heavy nucleus, so that the influ- 
ence of the potential attributed to the a particle needs only to be simulated. Of course, 
in this context, it is not assumed that the oscillator well is an appropriate tool to describe 
an alpha nucleus. So, to attain this purpose in a realistic manner, the superasymmetric 
two-center shell model (STCSM) described in Ref. || and improved in Ref. || was mod- 
ified in order to reproduce the single-particle levels assumed to describe an a cluster in 

the final stage of the disintegration. Consequently, the alpha oscillator stiffness hu 2a was 

1 /s 

forced to vary gradually from the usual STCSM value hu 2 = 41A 2 when the normalized 
coordinate R„—0, to the value extracted from Ref. [JIIJ huj a — 21.8 MeV when the scission 
point is reached, that means R n = 1. The nuclear shape parametrization being defined 
by two intersected spheres of different radii, Ri for the daughter and R 2 for the emitted 
fragment, the single generalized coordinate remains the elongation R, denoting the distance 
between the centers of the nuclei, or the normalized coordinate R n = (R — Ri)/(Rf — Ri), 
where Ri = Ro — R 2 , Rf = R% + Ri and Ro denotes the radius of the parent. Subsequently, 
another three modifications were introduced in the present version of the STCSM in order 
to obtain good a-level energies associated with the value huj a — 21.8 MeV. The first term 
characterizing the mass-asymmetry along the p-axis of Eq. (43) of Ref. 0] (proportional 
to hu 2 — Tuoi) was multiplied by the ratio (u 2a — uJi)/ {oj 2 — uJi). Another mass- asymmetry 
term along the z-axis was diagonalized using the eigenvalues of the asymmetric two-center 
potential: 
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In this way, energy shifts in single-particle levels associated to the change from u 2 to u)2 a 

1 /2 

are produced, so that a good two-center oscillator with fwoi = 41 A/ and fto^ = 21.8 MeV 
is obtained after the scission. Finally, the spin - orbit and L 2 interaction matrix elements 
associated to the light fragment were multiplied by the ratio uj 2a /u 2 in order to be coherent 
with the formal definitions of these interaction terms which assume proportionality to hu 2 . 
The signification of the notations can be found in Ref. 0. 

The spin-orbit and L 2 coefficients, together with the depths of the potentials for the 211 Po 
and 207 Pb used in this work are k = 5.75 x 10~ 2 , i] = 0.43 and ^=55.80 MeV, respectively. 
These values for the smooth spherical oscillator potentials were deduced from a fit of the 
208 Pb single-particle energies given in Ref. |T|: E lill/2 =-3.16 MeV, E 2gg/2 =-3.94 MeV, 



E 3pi/2 =-7.37 MeV, E 2h/2 =-7.94 MeV, E 3py2 =-8.36 MeV and E Ul3/2 =-9.00 MeV. The 



single-particle energies of 4 He are: -12 MeV, 1.38 MeV, 1.68 MeV, 5.34 MeV, 10.43 
MeV and 11.86 MeV for the levels lsi/2, 1^3/2, lpi/2, 2si/ 2 , ld 5 / 2 and 1g? 3 / 2 , respectively. A 
good fit for the lower energies in the frame of the oscillator model with hu a =21.8 MeV was 
found to be: k=4.51x10~ 3 , rj =0 and 14=52.91 MeV. 

The STCSM level scheme is plotted in Fig. |I[ In the following, the condition of consis- 
tency is also achieved, the same shape for the microscopical and phenomenological models 
involved. The parent and the daughter do not have pronounced deformations, so that their 
ground state nuclear shapes can be approximated with spheres. The ground state config- 



uration of 211 Po is ('x{hQ/ 2 ) 2 i>(g 9 / 2 ) 1 ) 9 / 2 + JT31. For clarity, the levels of the parent will be 
labeled with the superscript P, while the superscripts D and a will be used for the daughter 
and the a nuclei, respectively. Up to now, the model evidences the variation of the levels 
for a modification of the nuclear shape, and indicates the origin of the nucleons belonging to 
the a-particle. As an interesting feature, the STCSM predicts that the linked level ls^ 2 of 
the a particle emerges from the orbital lg$/ 2 of the spherical 211 Po, which is deeply located 
in the parent potential well. The present formalism intends to explain the fine structure by 
considering single-particle transitions due to the radial and the rotational couplings. The 
levels with the same good quantum numbers associated to some symmetries of the system 
cannot in general intersect, but exhibit quasi-crossings, or pseudo-crossings, or avoided level 
crossings. The system is characterized by an axial symmetry, therefore the good quantum 
numbers are the projections of the nucleon spin Q. The radial coupling causes transitions of 
the unpaired nucleon near the avoided level crossings. True crossings can also be obtained 
between levels characterized by different quantum numbers. Generally, the rotational cou- 
pling has a maximum strength in the vicinity of the true crossings. Transitions due to both 
couplings are taken into account in order to explain the excitations of the unpaired nucleon. 
It can be considered that the last unpaired neutron, initially located in the orbital 2g^ 2 
has the same chance to choose one of the levels with the projection Q included between 1/2 
and 9/2 when the nucleus starts to deform, therefore the occupation probabilities are the 
same. Comparing the diagrams (a)-(d) of Fig. 0, it is clear that only the level emerging 
from 2(^2 with the spin projection f2=l/2 finally reaches adiabatically the 3p^ 2 ground- 
state of the daughter, which means, a disintegration performed without excitations. Due 
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to the rotational coupling, a nucleon initially in the state 2g^ 2 ^=3/2 can jump in the 
state Q=l/2 during the disintegration, contributing in a smaller measure to obtain finally 
the daughter ground-state. Even by taking into account the rotational coupling, the other 
levels with Q >3/2 emerging from 2g^ 2 have a negligible contribution in the ground-state 
channel. Moreover, if the Coriolis coupling is not taken into account (using only the radial 
coupling), the levels with Q ^1/2 emerging from 2g^ 2 attain finally higher daughter levels 
with at least 3 MeV excess in energy (for example the orbital 2g^ 2 of the daughter). As the 
penetrability decreases dramatically with the height of the barrier, it becomes clear that the 
processes characterized by a nucleon emerging from 2gg, 2 Q >3/2 in the initial moment of 
the disintegration are unlikely. This discussion allows us to fix the initial conditions for our 
process: initially, the valence nucleon of the decaying system can be considered in the state 
2gg, 2 with Q=l/2 and 3/2. To show how it is possible that the unpaired nucleon arrives on 
allowed excited states of the daughter due only to the radial effect, an arrow in Fig. |2| (a) 
indicates some avoided level crossings between the levels which reach adiabatically the 3pP, 2 , 
2/5/2 an d 3p^/ 2 daughter orbitals. The transition probabilities are strongly enhanced in the 
avoided level crossing regions in accordance with the Landau-Zener effect. The transition 
probability between two adiabatic single-particle states strongly depends on the velocity 
of passage through the avoided crossing regions v tn n, and, implicitly on the tunneling time 
of the barrier. When the velocity increases, the transition probability is enhanced and the 
nucleon follows with a larger probability the so called diabatic state. Using the above de- 
scription, it is intended to reproduce the fine structure |13| exhibited by the a decay of 211 Po: 
98.9 % transitions to the l/2~ (level 3p° 2 ) ground state of the daughter, 0.55% transitions 
to the 5/2~ (level 2fj^ 2 ) first single-particle excited state and 0.54% transitions to the 3/2~ 
(level 3pg/ 2 ) second single-particle excited state. 

As briefly mentioned, both radial and rotational couplings caused by the relative mo- 



tion between the nascent fragments |14| , |15H are taken into account in order to calculate 
the occupation probabilities of several levels of the daughter by the unpaired nucleon. For 
simplicity, the effect due to the radial coupling is considered to be well reproduced by the 
Landau-Zener promotion mechanism p6|-|l~8| in the avoided crossing regions. Therefore, the 
fine structure of the process is strongly related to the dynamic characteristics of the system. 
The rotational or Coriolis coupling causes transitions between two levels for which the value 
of Q differs by one unit and it is proportional to the angular momentum operator of the 
single-particle j± matrix element. The STCSM provides the ingredients for calculating the 
single-particle transitions probabilities due to the Landau-Zener effect and to the Coriolis 
couplings: the interaction energies between the diabatic states in the avoided crossing 
regions, the diabatic level energies ef k (using spline interpolations) and the wave functions 
required to compute (flk\j±\^k =F 1) as described in Refs. [f[4l,|15| . The behavior of these 
ingredients are plotted in Fig. |[ The relative velocity between the nascent fragments can 
also be calculated quantum mechanically using a method similar to that of the variation of 
constants ||19|| , but in the following it will be considered as a fit parameter as in Refs. |5|-0. 
To obtain the final occupation probabilities of the daughter levels by the unpaired nucleon, 
a system of differential coupled equations must be solved: 
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with a\j lUm = /o(C - e) tm )dt/k, B R is the effective mass along the elongation R which 
was taken approximatively equal to the reduced mass of the system and 1=9/2 is the total 
spin of the system. The time dependence in the above equation can be removed by using 
the relations cf k = v tun dcf k / dR and R = v tnn t. The coefficients (c^ fe ) 2 give the occupation 
probabilities of the diabatic levels {j, £l k }- To solve this system, following the above discus- 
sion and inspecting the Fig. 0, it was considered that it is sufficient to choose the initial 
conditions so that the levels with Q=l/2 and 3/2 emerging from 2g^ 2 have the same initial 

occupation probabilities, which means that the equality (c^ 2 ) 2 + (^9/2) 2 = 1 is fulfilled. 
Also, by solving the system (|0.2j) , it is satisfactory to take into account the levels with £1=1/2 
emerging from lz p 1/2 , 2# P /2 , 3pf /2 , 2/ p /2 , 3p P /2 , those with £1=3/2 emerging from 2# P /2 , 2/ 5 p /2 , 
3p P / 2 , and those with £1=5/2 emerging from 2g^ 2 , fi p 3 / 2 . All the avoided crossing levels be- 
tween the selected adiabatic states are taken into account. Levels with £1 >5/2 do not reach 
the final channels we are interested in. 

For each channel, the penetrability PL im _(Qi) of the barrier was obtained using the nu- 
merical superasymmetric fission model |[20|| , the nuclear part being given by the Yukawa- 
plus-exponential approximation. This penetrability depends on the Qj-value of the channel 
i (i labels here the single-particle state of the daughter) and of the relative motion orbital 
momentum L im .. In the final channel 3pw 2 , due to the conservation laws, L im . has the value 
5 (rrii=l), in the final channel 2/° 2 , Li mi can be either 3, 5 or 7 (mj=l,2,3) and in the final 
channel 3pg/ 2 , L imi can be either 3 or 5 (mj=l,2). For a specific final single-particle state, 
for example 2/B 2 , it is not possible to discriminate between the possible values of the relative 
motion orbital momentum Lj m4 =3,5 and 7 in order to compute only one barrier penetrability. 
In these circumstances, by analogy to the Mang's formulae of Refs. 0|| concerning the ra- 
dial motion and the associated angular momentum, we consider that the angular momentum 
L im - , used in calculating the penetrabilities, has a probability to be obtained in the final chan- 
nel directly proportional to the square of the Clebsh-Gordon coefficient (jl£l — £l\L im .O) 2 . 
So that, the spectroscopic amplitude in the channel i associated to the spin £l k and the 

momentum L im! . will be p { "™ ,; = (c* k ) 2 (jil£l k - £l k \L im >0) 2 / T, mi {jil^k ~ £lk\L imi 0) 2 where 
the summation on m ; is done on the allowed values of L im .. The partial half-life T^ fc for 
the channel {i, £l k } becomes proportional to the quantity: 

T? k oc , (0.3) 

Em t Pi Ph im . (Qi) 

the proportionality factor being given by the barrier assault frequency. The partial half- 
lives for the transitions to the ground-state To„d , to the first excited state T 9f D and to the 

° Jp l/2' ^5/2 

second excited state TLd are: 

•^3/2 
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The barrier assault frequency being the same for all the channels, the relative intensities 
T 3p D ^/T 2 fD^ and T 3p D jT 3p o ^ for the fine structure can be obtained. Several tunneling 
velocities, considered here as a fit parameter, have been tried. For a tunneling velocity 
of 9xl0 6 fm/fs, the ratio between the intensity for transitions to the first excited state and 
to the ground state was found to be 0.0071 and the obtained ratio of the same parameter 
between the second excited state and the ground state was 0.0062. These results are in good 



agreement with the experimental values presented before. Moreover, calculations of Ref. fl9 
show that in the quantum time-dependent approach, the tunneling velocity is of the order 
of lxlO 7 fm/fs. These calculations suggest that the a decay fine structure phenomenon 
can be explained quantitatively by describing the decaying system with molecular models, 
and it can be stated that the quantitative characteristics of this phenomenon are ruled by 
dynamical effects. In an avoided crossing region, the two eigenf unctions of the adiabatic 
levels exchange their characteristics. If the relative distance R change infinitely slow, the 
unpaired nucleon will remain in the same adiabatic single particle state after the passage 
through a quasi-crossing, any other available single particle excited state being unfavoured. 
For a large tunneling velocity, the unpaired nucleon will follow the diabatic state after the 
passage through a quasi-crossing, all the other states being unfavoured. For a finite tunneling 
velocity, an intermediate situation arises, some single particle states being favoured and other 
single particle states being unfavoured. The model propose here, as the usual picture that 
consists of calculating the overlaps between the parent and the channel wave functions, is also 
based on the existence of favoured and unfavoured transitions in odd nuclei. The proposed 
formalism offers a competitive description of the alpha decay mechanism, by investigating 
for the first time the modality in which the levels initially bunched in shells are reorganized 
during the disintegration to realize the final energy configuration. 
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FIGURES 



FIG. 1. Neutron level scheme for the a decay of Po as function of the normalized elongation. 
The spectroscopic notations are used to describe the levels of the parent nucleus in the left side of 
the figure while in the right side of the figure the levels of the daughter and the alpha particle are 
labeled in the first and second column, respectively. 

FIG. 2. Detailed part of the level scheme. The levels with = 1/2 (a), SI = 3/2 (b), Si = 5/2 
(c), Si = 7/2 (d), Si = 9/2 (e) and Si = 13/2 (f) are plotted with thick lines. The levels with 
Si = 11/2 (e) and Si = 15/2 (f) are plotted with thick dotted lines. 

FIG. 3. Differences between the adiabatic levels with spin Sl=l/2 are presented in the plots 
(a) to (f) suggesting the possible avoided crossing regions, that means, points of nearest approach 
between two adiabatic levels. These regions are marked with arrows and the spectroscopic notations 
are displayed on each plot. The adiabatic Ei and the diabatic ej levels are presented in the picture 
(g) with full and dotted lines respectively only in the case Sl=l/2. The interaction energies 
between the diabatic levels of the plot (g) are presented in (h). The levels with spin Sl=\/2 
emerging from 2/f, 2 , 3p 3 / 2 and l^ 3 / 2 are presented as full lines while the level with spin 3/2 
emerging from li^ 3 / 2 is plotted as dotted line in picture (i). For the levels presented in (i), the 
matrix elements (3/2|j + |l/2) = (l/2|j_|3/2) are drawn in (j). Asymptotically, for R n =0, the 
matrix element between Sl=l/2 and 3/2 of the levels belonging to the subshell l^ 3 / 2 (f un line) 
has the value (SI ± l|j'±|0) = hy/Jj TSl)(j ±Sl + 1) while for R n — > oo the same value is obtained 
within the matrix element for the level Sl=l/2 emerging from 3p 3 ^ 2 and the level lif 3/ / 2 Sl=?>/2 
(dashed line). Otherwise, asymptotically the values are zero (dotted line). A pronounced maximum 
of the matrix elements is obtained when the levels with different spin projections intersect. 
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